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THE EFFICIENCY OF GRAIN ALIGNMENT IN DENSE 
INTERSTELLAR CLOUDS: A REASSESSMENT OF 
CONSTRAINTS FROM NEAR INFRARED POLARIZATION 

D. C. B. Whittet 1 , J. H. Hough 2 , A. Lazarian 3 & Thiem Hoang 3 

ABSTRACT 

We present the results of a detailed study of interstellar polarization effi- 
ciency (as measured by the ratio p\/r\) toward molecular clouds, with the aim 
of discriminating between grain alignment mechanisms in dense regions of the 
interstellar medium. The data set includes both continuum measurements in the 
K (2.2 jum) passband and values based on ice and silicate spectral features. Back- 
ground field stars are used to probe polarization efficiency in quiescent regions 
of dark clouds, yielding a dependence on visual extinction well-represented by a 
power law (jp\ j T\ oc [Ay]~ ' 52 ), in agreement with previous work. No significant 
change in this behavior is observed in the transition region between the diffuse 
outer layers and dense inner regions of clouds, where icy mantles are formed, 
and we conclude that mantle formation has little or no effect on the efficiency 
of grain alignment. The field-star data are used as a template for comparison 
with results for embedded young stellar objects (YSOs). The latter generally 
exhibit greater polarization efficiency compared with field stars at comparable 
extinctions, some displaying enhancements in p\jr\ by factors of up to ~ 6 with 
respect to the power- law fit. Of the proposed alignment mechanisms, that based 
on radiative torques appears best able to explain the data. The attenuated ex- 
ternal radiation field appears adequate to account for the observed polarization 
in quiescent regions for extinctions up to Ay ~ 10 mag. Radiation from the 
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embedded stars themselves may enhance alignment in the lines of sight to YSOs. 
Enhancements in p\/r\ observed in the ice features toward several YSOs are of 
greatest significance, as they demonstrate efficient alignment in cold molecular 
clouds associated with star formation. We advocate spectropolarimetry of ice 
features in background stars as a reliable means of mapping internal magnetic 
fields within molecular clouds. To date, such observations exist for just one line 
of sight. 

Subject headings: polarization — dust, extinction — ISM: magnetic fields - 
infrared: ISM 



1. Introduction 

The detailed physics of grain alignment is yet to be fully understood despite much activ- 
ity and considerable progress in the years since the pioneering work of Davis & Greenstein 
(1951; DG) (see reviews by Lazarian 2003, 2007). Developing a comprehensive, realistic 
model for grain aligment is an important astrophysical goal for several reasons. Magnetic 
fields play vital roles in interstellar processes and star formation, and the polarization that 
results from grain alignment provides the means to study magnetic fields if the alignment 
process is well understood. Polarization is also a valuable tool for studying the dust grains 
themselves, placing constraints on their magnetic properties, structure, shape and size distri- 
bution. It is generally assumed that the dominant process responsible for alignment involves 
interaction between the spin of the particles and the ambient magnetic field, as originally 
proposed by DG: a spinning grain tends to become orientated with its longest axis perpen- 
dicular to the angular momentum vector, and paramagnetic relaxation imposes alignment 
of the angular momentum with respect to the magnetic field vector, B. For polarization 
produced by dichroic extinction, the mean direction of the electric vector in the transmitted 
beam is thus parallel to B, and so, in principle, the observed polarization traces the magnetic 
field in the plane of the sky. In practise, the observed polarization depends on the mean 
field direction, weighted according to the distribution of dust along the line of sight and the 
efficiency with which it is aligned. 

The original DG formulation assumed the grains to be spinning thermally and to have 
paramagnetic properties (see Roberge & Lazarian 1999 for a modern account). This mecha- 
nism was qualitatively successful, explaining the pattern of interstellar polarization observed 
in the Milky Way, but it failed quantitatively to predict significant alignment under realistic 
physical conditions. Subsequent investigations considered enhancements to the DG mecha- 
nism that might arise if the grains are endowed with superparamagnetic properties or if they 
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spin suprathermally. The Purcell (1979) mechanism attributes suprathermal spin primarily 
to H 2 formation at active sites on the grain surfaces. More recent work has emphasized ra- 
diative torques (RT; Dolginov & Mytrophanov 1976) as a promising means of driving grain 
alignment (Draine & Weingartner 1996, 1997; Weingartner & Draine 2003; Cho & Lazarian 
2005). An analytical model of alignment by RT proposed by Lazarian & Hoang (2007) gives 
physical insight into the nature of the mechanism, showing that alignment occurs with long 
grain axes perpendicular to B, i.e. in the expected orientation; this results not from para- 
magnetic relaxation, however, but from averaging of RT as the grain angular momentum 
precesses about B. 

An important objective from the observational standpoint is to attempt discrimination 
between these proposals by studying interstellar polarization in environments where one or 
other mechanism may be expected to fail. Key questions include how the efficiency of align- 
ment varies with gas density and with the intensity and spectral energy distribution of the 
ambient radiation field. The polarization of starlight at visible wavelengths provides infor- 
mation on alignment in predominantly diffuse, translucent phases of the ISM. Observations 
in the near infrared allow stars behind much larger columns of dust to be explored, including 
those along lines of sight that intercept dense molecular clouds. Results have shown that the 
polarization per unit extinction declines systematically with optical depth in such regions 
(Goodman et al. 1992, 1995; Gerakines et al. 1995), implying that grains within clouds are 
generally less well aligned than those in the diffuse (intercloud) ISM. This behavior is also 
indicated by observations of polarized far infrared and submillimeter emission, which show 
evidence for reduced (but non-zero!) polarization amplitudes toward the centers of isolated 
globules and starless core (Hildebrand et al. 1999; Ward-Thompson et al. 2000; Wolf et al. 
2003). At first sight, this behavior does not appear to allow effective discrimination between 
alignment mechanisms, as it is qualitatively consistent with all of them (Lazarian et al. 
1997). In the absence of embedded sources, alignment by radiative torques will naturally 
decline in efficiency within a dark cloud as the external radiation field becomes attenuated. 
Spin-up by Purcell torques will decline as H — > H2 conversion proceeds toward completion at 
higher densities, and is moreover rendered ineffective by resurfacing of grains (as evidenced 
by the appearance of ice mantles) which neutralizes the active centers at which H recombi- 
nation preferentially occurs (Lazarian 1995)B Finally, paramagnetic or superparamagnetic 



lr This might be explained, at least in part, by the relatively low spatial resolution of the far infrared 
observations if changes in magnetic field direction occur in the cores on scales smaller than the beam size. 

Subsequent work by Lazarian & Draine (1999a, b) showed that grains with sizes smaller than ~ 1 fjxa 
flip rapidly, averaging out Purcell's torques. This effect, termed "thermal trapping", diminishes the effects 
of uncompensated torques arising from H2 formation. 
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alignment at thermal rotation speeds will fail as the dust and gas temperatures equilibriate 
at high density. Indeed, the question can be turned around: which mechanism, if any, is 
capable of sustaining the observed (albeit reduced) levels of grain alignment within dense 
molecular clouds and starless cores? 

The puzzle of grain alignment within starless cores was addressed by Cho & Lazarian 
(2005), who noted that RT alignment is most effective for the largest grains (mean radii 
> 0.3 /im) in the size distribution, which are aligned efficiently by near infrared radiation 
that penetrates much deeper into a cloud than ultraviolet or visible light (see also Bethell 
et al. 2007). The RT mechanism thus predicts alignment of the grains responsible for near 
infrared extinction and polarization at substantial optical depths within dense clouds, where 
other mechanisms are predicted to fail. 

Given the importance of magnetic fields and the potential of far infrared observations for 
mapping them in dense phases of the ISM (e.g. Wolf et al. 2003; Bethell et al. 2007; Pelkonen 
et al. 2007), it is timely to reexamine observational constraints placed by the polarization of 
starlight at shorter infrared wavelengths. The purpose of this paper is to examine in detail 
the correlation of alignment efficiency with extinction through representative dense clouds, 
with the goals of (i) discriminating between alignment mechanisms, and (ii) reassessing the 
value of near infrared polarimetry as a tracer of the internal magnetic fields. The polarization 
efficiency, i.e. the ratio of linear polarization amplitude to optical depth, p\/r\ (also referred 
to as the specific polarization), is taken as an observational measure of alignment efficiency, 
subject to caveats described in §2.1 below. This approach is not new (cf. Gerakines et al. 
1995) but has the advantage that the Two-Micron All-Sky Survey (2MASS) is now available 
to provide a coherent and essentially complete photometric database for making extinction 
estimates toward stars in the polarization database. We also utilize polarization data pub- 
lished since the earlier study, including observations of ice and silicate features as well as the 
near infrared continuum. 



2. Data 
2.1. Rationale 

The value of polarization efficiency p\/r\ observed at some wavelength A toward a 
background star depends on several factors in addition to the desired one, i.e. the efficiency 
with which the grains along the line of sight are aligned by the ambient magnetic field. For 
perfect alignment in a uniform field B, P\/t\ changes with the viewing angle relative to 
B (no polarization is seen when the line of sight is parallel to the magnetic field). If field 
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lines are twisted, or if two or more discrete clouds or clumps are present in the line of sight 
with different field orientations, then the net polarization is reduced relative to extinction. 
These unwanted effects can be minimized by careful target selection: well-sampled, isolated 
clouds with negligible foreground or background extinction and permeated by a uniform 
magnetic field with a strong transverse component are ideal. The well-known system of 
dark clouds in Taurus is the best available approximation: it is nearby (140 pc), lies in a 
direction of little general extinction, and is reasonably well-sampled by background field 
stars exhibiting substantial degrees of interstellar polarization (e.g. Whittet et al. 2001 and 
references therein). Moneti et al. (1984) used optical polarimetry of background stars to 
show that the macroscopic B field is, indeed, reasonably uniform across the cloud. 



2.2. Taurus 

Table 1 presents a compilation of broadband K (2.2 jttm) linear polarization and ex- 
tinction data for stars subject to reddening by dust in the Taurus dark cloud. Most are 
background field stars that lack any direct association with the cloud; for comparison pur- 
poses, available data for a few embedded young stellar objects (YSOs) that meet our selection 
criteria (below) are also included in Table 1. Additional spectropolarimetric data and optical 
depths available for ice features (H 2 and CO) in one Taurus field star (Elias 16) are listed 
in Table 2. In general, polarization data were accepted only if the fractional error in the 
amplitude (Ap/p) is less than 25%. Exceptions to this quality filter were made in a few cases 
(identified in the footnotes in Table 1) where polarimetry at shorter wavelengths (e.g. the J 
and H bands) could be used to corroborate the K value using the usual empirical form for 
the wavelength-dependence of polarization (e.g. Whittet et al. 1992 and references therein). 
In all cases where multiple-passband polarimetry is available, a check was made to ensure 
that the spectral dependence is broadly consistent with interstellar polarization dominated 
by a single cloud. 

Infrared photometry for all selected objects was collected from the on-line 2MASS cat- 
alog and is generally accurate to 0.02 mag or better. Two methods were used to evaluate 
visual extinction (Ay) and K-band optical depth (tk), dependent on the additional data 
available in each case. For objects with reliable visual photometry and spectral classifica- 
tion^ (typically those with Ay < 5 mag), Ay is estimated using the relation 

Ay « I.IEv-k (1) 



3 Spectral classifications are available from the SIMBAD database and references therein for all stars in 
Table 1 with Elias or HD prefixes, and also for Tamura 8. 
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where E V -k is the V — K color excesses arising from interstellar reddening. The factor 1.1 
in eq. (1) effectively represents extrapolation of the mean interstellar extinction curve from 
2.2 /im to the long wavelength limit (Whittet 2003), and appears to be a robust value not 
strongly dependent on the assumed grain model over the usual range of observed extinction 
curves (Whittet & van Breda 1978). This u Ey-K method" (Table 1) provides Ay values 
accurate to ±0.1 or better: most of the results that use it are taken from Whittet et al. (2001), 
who provide further discussion and references to spectral classifications, visual photometry 
and intrinsic colors. The fT-band optical depth then follows from the relation 

t k = Ajr/1.086 = (A v - E V _ K )/1.086 (2) 

which in combination with eq. (1) yields t k « 0.084Ay. 

For objects in our sample lacking visual photometry, we rely on infrared photometry 
alone to evaluate 

A v = ri Ej_ K (3) 

and 

r K = r 2 Ej_ K (4) 

(the u Ej_k method"), where r\ and r 2 are factors that depend on the form of the extinction 
curve over the relevant wavelengths (Martin & Whittet 1990; Whittet et al. 2001). Mean 
values r\ = 5.3 ± 0.3 and r 2 = 0.50 ± 0.03 determined from data for the Taurus cloud 
(Whittet et al. 2001) are adopted here (for comparison, r\ 6.0 and r 2 ~ 0.54 in the diffuse 
ISM). For stars with known spectral classifications, Ej_ K is determined routinely from the 
observed colors to yield Ay values accurate to ±0.4 or better. For stars lacking spectral 
classifications (A v values followed by colons in Table 1), Ej_ K = Ej_ H + E H _ K is estimated 
from the observed locus in the J — H vs. H — K diagram (Fig. 1) by extrapolation along the 
appropriate reddening vector onto intrinsic color lines. This generally provides unambiguous 
results because background stars bright enough to be included in our sample are expected 
to be predominantly either late-type (K, M) giants or main-sequence stars of earlier type 
(typically B-G); red dwarfs distant enough to be background to the cloud are too dim at 
2.2 fxm for polarimetry to be feasible with current technology. Obviously, the absence of 
spectral type information reduces the reliability of Ay estimates (±1 mag would be a worst 
case scenario) but the fractional error is nevertheless quite low for stars with the highest 
extinctions (Ay > 10 mag). Note also that this technique improves upon previous estimates 
that adopted a standard "one fits all" intrinsic J — K color for background stars lacking 
spectral classifications (Tamura et al. 1987; Goodman et al. 1992; Gerakines et al. 1995). 

The distribution of Taurus YSOs in the J — H, H — K color-color diagram is compared 
with that of the field stars in Fig. 1. As expected, the field stars occupy a zone consistent 



- 7- 



with normal intrinsic colors, subject to interstellar reddening. The trend for YSOs is similar, 
but subject to a characteristic displacement to the right that signals a contribution to their 
infrared colors from circumstellar matter (e.g. Itoh et al. 1996). Our estimates of Ay for 
these objects assume that the infrared excess affects predominantly the H — K color and is 
negligible in J — H. 

2.3. Ophiuchus 

Although not as well sampled as Taurus for near infrared polarimetry, the p Oph Dark 
Cloud provides a valid and interesting comparison. Both clouds are regions of low-mass 
star formation, but p Oph is subjected to a harsher radiative environment arising from its 
proximity to an OB association (Sco OB2): Liseau et al. (1999) deduce that the local UV 
field in p Oph is typically a factor of 20-140 above the interstellar average. 

Table [3] presents a compilation of available polarization data, together with extinction 
estimates, for stars reddened by dust in this region, using the same selection criteria as 
before (§2.2). Both background field stars and embedded YSOs are included. The color 
excess methods used to estimate Ay are the same as those described in §2.2, except that 
values of the parameters T\ = 4.9 ± 0.3 and r 2 = 0.53 ± 0.05 (eq. 3 and 4) measured from 
the mean p Oph extinction curve (Martin & Whittet 1990; Kenyon et al. 1998) are adopted 
in place of the Taurus values. The color-excess method is unreliable for deeply embedded 
objects such as the Class I YSO Oph 29 because its near infrared colors are dominated by 
circumstellar emission; in this case, Ay is estimated from the depth of the observed 9.7 /im 
silicate absorption: 

A v « 18r 9 . 7 (5) 

(Whittet 2003). This method assumes that the observed correlation between t$j and Ay in 
the diffuse ISM holds inside the dense cloud. Note, however, that a recent study by Chiar 
et al. (2007) suggests a systematic difference in the sense that eq. 5 may underestimate the 
true Ay. 

2.4. Other datasets 

Several additional sets of polarimetric data that provide useful comparisons for those 
discussed in detail above are listed in Table HI For the most part, these are self-explanatory 
and accessible from the single reference listed; data for individual objects are not included 
here. LI 755 is a filamentary dark cloud cloud adjacent to and physically associated with the 
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p Oph cloud (Loren 1989), and L1641 is an active site of star formation in the Orion region: 
all point sources listed in the cited references that meet our criteria were selected. The other 
entries in Table @] are grouped by spectral feature (H 2 0-ice, CO-ice, silicate) rather than by 
region of sky. For a given absorption feature, the same observation provides information on 
the peak values of p\ (from Stokes parameters Q and U) and t\ (from intensity I), measured 
in each case relative to the continuum, leading to a measure of p\/r\ specific to the particles 
responsible for the feature (see Aitken 1996 for a review). This is straight-forward for ice 
features, which are only ever observed in absorption. However, silicate features may appear 
in absorption, emission, or a superposition of the two, and care was taken to select only 
those in which the observed 9.7/zm profile is dominated by absorption (Smith et al. 2000). 
Table [5] presents a comparison of polarization efficiencies for YSOs determined from ice and 
silicate features and i^-band continuum data. 

The methods used to estimate Ay are exactly analogous to those discussed above. In 
the case of L1755, the p Oph extinction parameters (§2.3) are assumed. In the case of L1641, 
we adopt values based on H — K colors from Casali (1995). Finally, Ay values calculated 
from the silicate optical depth use T9.7 values from Smith et al. (2000). 

3. Results and Discussion 

3.1. The Taurus Cloud 

The analytical tool used in this study is the plot of polarization efficiency vs. visual 
extinction (p\/r\ vs. Ay] Figs. 2-8). For lines of sight passing through an idealized, isolated 
cloud consisting of denser clumps within a more diffuse envelope, Ay is expected to be a 
reasonable proxy for mean density because stars that lie behind clumps will have systemat- 
ically higher extinctions. This assumption is corroborated by the fact that ices are detected 
only in lines of sight with Ay above some threshold value (A^ « 3.2 mag in the case of the 
Taurus cloud; Whittet et al. 1988, 2001). As ices do not survive in diffuse regions regions 
but appear to be widespread inside molecular clouds, this extinction threshold is assumed 
to correspond approximately to the transition between atomic and molecular phases of the 
ISM (Whittet et al. 2001). 

Plots of P\/t\ vs. Ay for Taurus stars are shown in Fig. 2 (see Table [6] for a key to 
plotting symbols in all figures). The complete set of K-band data for field stars (Table 1) 
is included in frame (a), and the least-squares power-law fit to this dataset (see below) is 
shown in all frames (thick diagonal line). The thinner (parallel) diagonal lines, which lie 
0.34 dex above and below the best fit, represent the approximate upper and lower bounds of 
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the distribution. Frames (b)-(d) in Fig. 2 show subsets of the data in frame (a) according 
to the method used to determine Ay and r K (E V _ K method; Ej_ K method without or 
with spectral types, respectively). The Ay range spanned by these subsets is different, as 
expected, but they nevertheless show reasonable consistency with the common fit, confirming 
that no major systematic errors exist between them. It is also clear that there is real scatter 
in the data: the degree of scatter is significantly greater than that implied by observational 
error, at least for the stars with spectral classifications. This is to be expected for reasons 
explained in §2.1. 

Fig. 2a shows that a power law provides a fair representation of the correlation between 
Pk/tr and Ay for field stars, in agreement with the results of Gerakines et al. (1995). The 
fit is of the form 

Pk /t k = a(A v )- b (6) 

where a = 5.3±0.3 and b = 0.52±0.07. Power-law behavior is consistent with the dependence 
of magnetic field strength on density for classical DG alignment (Vrba et al. 1981), although 
the applicability of such a model is questionable (see Gerakines et al. 1995 for discussion 
and caveats). We use it here merely as an empirical law that can be used to predict (by 
extrapolation) the alignment to higher extinctions. 

Fig. 2e compares the distribution of 2.2 /im broadband data for Taurus YSOs (Table 1) 
with the fit to field stars. A systematic difference would be expected if, for example, radiation 
from embedded stars influences grain alignment in surrounding regions of the cloud. There 
is, indeed, some evidence for higher polarization efficiency toward the YSOs (6 out of 7 lie 
above the correlation line) but the sample is too small to make a strong statement. 

Fig. 2f plots two Px/tx values based on ice-feature spectropolarimetry of Elias 16 (Ta- 
ble 2). Polarization is measured at the peak of each feature relative to the continuum, and 
ratioed with the corresponding optical depth (Hough et al. 1988, 2008). The results are 
consistent with the K-b&nd continuum datum for this star (Table 2) and with the general 
trend implied by the fit to all field star data in Fig. 2. This is an important result as it 
implies that ice-mantled grains in the line of sight are aligned to a similar degree compared 
with those responsible for the continuum polarization, even though the latter must include 
a contribution from unmantled grains. The detection of polarization in the CO ice feature 
is an especially significant result (Hough et al. 2008; Crysostomou et al. 1996), as the low 
sublimation temperature of CO (~ 17 K in pure form) implies that it can exist in solid phase 
only in the coldest, densest regions along the line of sight. 

The vertical dotted line in each frame of Fig. 2 denotes the threshold visual extinction 
(A th = 3.2 ±0.1 mag) for detection of H 2 0-ice in the Taurus cloud, as discussed above. Below 
this threshold the grains resemble those in the diffuse ISM, lacking a volatile component; 
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above it they acquire icy mantles, as evidenced by the onset of the 3.0 /zm H 2 0-ice absorption 
feature in the spectra of background stars. CO2 and CO ices are also observed with somewhat 
higher thresholds (Ah ~ 4.3 and 6.7 mag, respectively; Whittet et al. 2007; Chiar et al. 
1995). A key question is whether the ability of a grain to align in the ambient magnetic 
field undergoes any substantive change with the growth of a mantle. The distribution of 
data (e.g. Fig. 3) shows no evidence for a significant discontinuity or change in slope near 
the threshold value of Ay. This general result, together with the presence of polarized ice 
features in at least one line of sight through the cloud, clearly demonstrates that grains 
are, indeed, capable of being aligned in environments where they are subject to ice-mantle 
growth. As the accretion process that leads to mantle formation is expected to inactivate 
the sites on a grain surface that drive Purcell alignment by H recombination (§1), this 
finding constitutes further evidence against the Purcell mechanism as the dominant means 
of alignment for extinctions Ay > Ah- 

The general trend of declining polarization efficiency with extinction led some investi- 
gators to propose that grains within dense clumps make no contribution to the net observed 
polarization of background stars (Goodman et al. 1992, 1995; Arce et al. 1998). Obser- 
vations of ice-feature polarization provide counter-examples but it is nevertheless useful to 
consider this scenario as a limiting case, which we refer to as the "skin-depth" model. The 
structure of the Taurus cloud system resembles a diffuse screen with embedded clumps, as 
described above, such that the maximum opacity of the screen is approximately equivalent 
to the ice threshold extinction (see, e.g., Fig. 1 of Whittet et al. 2004). Suppose that all 
of the observed .fT-band polarization arises in the diffuse screen at extinctions in the range 
< Ay < Ah (Ah = 3.2 mag). Then, for lines of sight that intersect clumps (Ay > Ah) 
Pk is independent of Ay, and thus Pk/tk oc [^4y] _1 (assuming tk oc Ay). This form of 
behavior is compared with the data in Fig. 3, the constant of proportionality being set by 
the mean of all data below the ice extinction threshold. A power- law index of —1 obviously 
represents a much steeper decline than our best fit (—0.52 ± 0.07). The skin-depth model is 
nevertheless reasonably consistent with the data for Ay < 10 mag (Fig. 3), failing at higher 
extinction: at Ay ~ 25 mag, for example, it predicts Pk/^k to be a factor of ~ 3-4 less than 
observed. 

The "Taurus trend" represented by the best fit and upper and lower bounds in Fig. 2 is 
used as a template for comparison with the other samples discussed in the following sections. 
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3.2. The Ophiuchus Cloud 

Polarization efficiency vs. visual extinction is plotted for p Oph stars in Fig. 4. Broad- 
band K data are from Table 3, and spectropolarimetric data for YSO Oph 29 (H^O-ice and 
silicate) are from references listed in Table 4. The diagonal lines show the fit and upper/lower 
bounds to Taurus stars, as in Fig. 2. A fit to p Oph field stars (not shown) yields a power-law 
index —0.45 ± 0.10, not significantly different from the Taurus result. In general, both field 
stars and YSOs in p Oph are consistent with the trend implied by the Taurus data. No 
systematic difference is apparent between p Oph field stars and YSOs in the Ay-range of 
overlap (~ 5-15 mag) in Fig. 4. At higher extinction {Ay ~ 15-30 mag), there is evidence for 
greater polarization efficiency in the YSO data. Interestingly, the H 2 0-ice datum (triangle) 
for Oph 29 is significantly lower (and consistent with the Taurus line) compared with the 
silicate and i^-continuum data in the same line of sight (see also Table [5]). 

3.3. L1755 

Polarization efficiency vs. visual extinction for L1755 is plotted in Fig. 5. All points 
represent field stars behind this filamentary dark cloud (Goodman et al. 1995). The diagonal 
lines show the fit and upper/lower bounds of the Taurus data, as before. The decline in 
polarization efficiency with increasing extinction appears to be significantly steeper in LI 755 
compared with Taurus (or Ophiuchus), and is more consistent with the skin-depth scenario 
described in §3.1 above, illustrated by the dotted line in Fig. 5. This finding is in agreement 
with the results of Goodman et al. (1995). 

It is not clear why L1755 should display this form of behavior, which is unique amongst 
the regions we investigated. It might be that physical conditions leading to grain alignment 
are lacking inside this cloud, that the internal grains are effectively spherical, or that pe- 
culiarities in magnetic field structure and/or viewing geometry conspire to systematically 
reduce the observed j>\/t\. It should also be noted that the total range in Ay covered by the 
observations (2-8 mag) is small compared with that of the other clouds studied: it would be 
of obvious interest to obtain polarization data for stars with Ay > 10 mag toward L1755 to 
see if the trend continues. 



3.4. L1641 

L1641 is a turbulent region of clustered low-mass star formation (e.g. Stanke & Williams 
2007), and as such it differs markedly from the prototypical dark cloud. Sources with avail- 
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able polarimetry (Casali 1995) are embedded YSOs originally detected by IRAS (Strom et al. 
1989). The level of polarization efficiency observed in L1641 is nevertheless remarkably con- 
sistent with that in our template cloud, as illustrated in Fig. 6. All 15 objects in the sample 
lie within the bounds implied by the Taurus data. 

3.5. Feature polarization in YSOs 

Young stellar objects with feature polarization data available are listed in Table O Out 
of 19 in the set, all but two (Oph 29 and SVS 13) are considered to be high-mass objects that 
seem likely to exert a strong influence on physical conditions in their vicinity. Fig. 7 plots 
feature polarization data for ices and silicates in frames (a) and (b), respectively, together 
with available -fT-band data for the same objects in frame (c). The plot for ices includes one 
datum for CO, the other points being for H2O. In each frame, the locus of each individual 
point may be described as either "normal" or "enhanced" in terms of polarization efficiency 
with respect to the Taurus template distribution, "enhanced" meaning significantly above 
the upper bound. YSOs with enhanced polarization efficiencies according to this definition 
are identified by values in bold type in Table [51 Some objects, notably Mon R2 no. 2 and 
OMC-1 BN, show similar behavior at all wavelengths. Others show discrepancies that seem 
likely to reflect differing physical conditions in the individual lines of sight, perhaps leading 
to differences in the distributions of volatile and non-volatile grain material as well as in the 
efficiency of alignment. 

The ice features are more likely to reflect conditions within the molecular clouds, whereas 
the silicate features are more likely to include a circumstellar component. The fact that 
enhancements are observed for ices in 6 out of 9 lines of sight toward YSOs (and 5 out of 
7 toward high-mass YSOs) therefore seems highly significant: this result implies that grain 
alignment is more efficient in molecular clouds that contain regions of active star formation 
(especially massive active star formation), compared with those that do not. Independent 
support for this conclusion is provided by observations of continuum circular polarization, 
the highest amplitudes of which are likewise found in lines of sight to massive YSOs (e.g. 
Buschermohle et al. 2005, Clayton et al. 2005, and references therein), and are best explained 
by models that assume efficient grain alignment in the envelope of the YSO (Whitney & Wolff 
2002). 



-13- 



3.6. Overview 

Fig. 8 displays summary plots that include data for all regions, passbands and features 
in this study, with field stars and YSOs plotted in frames (a) and (b), respectively. The field- 
star data (Fig. 8a) are consistent with a power-law dependence for polarization efficiency on 
extinction in the observed range < Ay < 25, with scatter of approximately a factor 2 about 
the mean relationship for the template cloud. Results for field stars provide unequivocal 
evidence that ice-mantled grains are aligned inside dark clouds, and that the onset of mantle 
growth has no dramatic effect on net alignment. 

The YSO data (Fig. 8b) are distributed either close to or significantly above the mean 
power-law relationship for field stars behind the template cloud. YSOs typically lie above 
field stars at comparable extinctions within the range of overlap (1 < Ay < 25), and en- 
hancements in P\/t\ by factors of up to ~ 6 with respect to the power-law fit occur toward 
some YSOs (predominantly those of high mass and with extinction Ay > 20 mag). En- 
hancements of this order are seen in the ice features as well as in the silicate feature and 
.fT-band continuum. The simplest explanation for the distribution of YSO data is that each 
line of sight contains a "quiescent cloud" component that follows the general trend, and 
a "modified" component in which alignment is enhanced by local conditions, the relative 
contributions of each differing from one line of sight to another. It is logical to assume that 
alignment in the modified component may be driven by radiation from the embedded stars 
themselves. 



4. Modeling 

We noted in §1 that a systematic decline in alignment efficiency with increasing density 
inside cold (starless) molecular clouds is a prediction of all the major alignment mechanisms 
(Lazarian et al. 1997). In regions cold enough for solid CO to accumulate on the grains 
(Chrysostomou et al. 1996) we expect catastrophic failure of paramagnetic or superparam- 
agnetic DG alignment at thermal rotation speeds. Significant alignment at suprathermal 
rotation speeds by the Purcell mechanism also appears to be ruled out (§1); however, align- 
ment by radiative torques seems more promising as a mechanism that can operate within 
molecular clouds, albeit with reduced efficiency as the external radiation field becomes pro- 
gressively attenuated with increasing Ay. Extinction by dust naturally changes the spectrum 
as well as the intensity of the ambient radiation field, effectively blocking the shorter wave- 
lengths whilst allowing the infrared to penetrate much further into a cloud. A key point that 
arises from theory (Cho & Lazarian 2005; Lazarian & Hoang 2007) is that the efficiency of 
the RT mechanism increases systematically with grain size and is optimized for grains with 
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dimensions comparable with the wavelength. Thus, the largest grains in the size distribution 
may be aligned by near infrared radiation in regions of a cloud that are optically thick to 
shorter wavelengths. 

The alignment of grains by radiative torques (Dolginov & Mytrophanov 1976; Draine & 
Weingartner 1996, 1997; Lazarian & Hoang 2007) is used to calculate the expected relation- 
ship between polarization efficiency and visual extinction through an idealized homogeneous 
cloud. We assume the cloud to have gas density = 10 4 cm" 3 , gas and dust tempera- 
tures T gas = 20 K and T^ ust = 15 K, radius r c i OU( j = 0.5 pc and extinction Ay = 12 mag 
(center to edge). Oblate spheroidal grains with axial ratio a/b = 0.5 and sizes following 
the MRN distribution (Mathis, Rumple & Nordseick 1977) are adopted with limiting sizes 
a min = 0.005 jum and a max = 0.3 /im. In calculating radiative torques, we take into account 
that real grains are irregular and can have significant helicity. Although crucial to RT, how- 
ever, irregularity and helicity are higher order effects when we deal with grain extinction. 
Indeed, it is customary to approximate irregular but swiftly rotating grains by spheroids, 
and in our simulations of light propagation we follow this custom. 

Radiative torque efficiencies were first calculated for irregular grains using the publicly 
available code DDSCAT (Draine & Flatau 1994) for the entire spectrum of the mean inter- 
stellar radiation field (0. 1-100 /xm; Mathis, Mezger & Panagia 1983). Grains of shape 1 in 
Draine & Weingartner (1996, 1997) were adopted, each consisting of 13 identical cubes, and 
the calculations were carried out over the full range of grain sizes. We then calculated the 
angular velocity of grains spun up by radiative torques arising from the attenuated radiation 
field for grain size in the range a min to a max , assuming perfect coupling of the grain axis of 
major inertia and angular momentum, and with both axes parallel to the radiation beam. 
Other alignment mechanisms such as DG and Purcell alignment (§1) were disregarded. The 
minimum critical size for alignment, a crit was then deduced for the criterion of suprathermal 
rotation, u; ra d/^thcrm > 3. This criterion is based on calculations in Hoang & Lazarian (2007), 
who studied the influence of gas bombardment on alignment by radiative torques using the 
Langevin equation approach: it was shown that when grains are aligned by radiative torques 
at an attractor point with an angular momentum greater than a value about 3J t h e rm, they 
are not randomized by gas bombardment. The polarization arising from aligned grains is 
given by 



where C is a normalization constant of the size distribution, R is the Rayleigh reduction 
factor that characterizes alignment (Roberge & Lazarian 1999), Q po \ — Q± — Q\\ is the 
polarization efficiency, and Q±,Q\\ are extinction efficiencies for the axis of major inertia 
perpendicular and parallel to the electric field, respectively. The corresponding optical depth 





is given by 



Tx = 



/ 



ds 




Ca 3 ' 5 ira 2 Q ext (a, A) da 



(8) 



where Q ext = [2Q± + Q\\)/3 is the total extinction. The polarization specified by eq. (7) is 
optimized when the magnetic field vector B is transverse to the line of sight, but in general 
the direction of B is arbitrary; this may be taken into account by replacing R in eq. (7) with 
R* = Rcos 2 (, where ( is the angle between B and the plane of the sky (see e.g. Fig. 3 of 
Roberge & Lazarian 1999). R* is thus a composite quantity dependent on the direction of B 
as well as on the degree of grain alignment characterized by the Rayleigh reduction factor. 

Systematic changes 111 G?<crit cl function of the radiation field leads to a prediction of 
the wavelength of maximum polarization (A max ) with respect to extinction: this prediction 
is shown and compared with observational data for the Taurus cloud in Fig. 9. As the 
interstellar polarization curve typically peaks at visible wavelengths, empirical determination 
of A max (from the Serkowski formula) requires availability of polarimetric data in the visible 
as well as the near infrared, and is thus possible only for lines of sight with relatively low 
extinction (Ay < 6 mag with one exception; Whittet et al. 2001). Considerable scatter is 
evident in Fig. 9, but the model reproduces the general trend in the data toward longer A max 
at larger Ay. Note that no changes in the size distribution with Ay are assumed in the 
model: variations in A max are a consequence of variations in the average size of the aligned 
component of the grains. 

The dependence of polarization efficiency on extinction predicted by the model is shown 
for several values of R* and compared with observational data for the Taurus cloud in Fig. 10. 
In addition to the near infrared ratio Pk/tr (Fig. 10a) we also plot p max /Ay (Fig. 10b). The 
advantage of the latter quantity is that the formal errors in the observational data are smaller 
(because p max is determined from a fit to several passbands instead of from just one) at the 
cost of a much smaller sample. Our calculations illustrate that the RT model is capable of 
producing alignment sufficient to explain the observed polarization efficiencies in the Taurus 
cloud for extinctions up to at least Ay = 10 mag. Indeed, the predicted rate of decline in 
p/r with Ay is less steep than the observed trend, especially in the near infrared (Fig. 10a). 
This may easily be explained by physical differences in the real cloud compared with the 
idealized model cloud, such as twists in magnetic field direction and superposition of discrete 
cloud components with different field orientations (§2.1). 

In our simple model, the alignment mechanism becomes inoperative as Ay approaches 
10 mag because even the largest grains fail to align (a crit — > a m ax)- However, RT alignment 
is still possible at higher extinctions if the grains are subject to growth (leading to an 
increase in a max ). The grains are, indeed, expected to grow inside the cloud by processes 
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of both mantle accretion and coagulation (e.g. Draine 1985; Wurm & Schnaiter 2002). It is 
important to realize, also, that Ay is not necessarily a reliable measure of radiation field, as 
it accounts for attenuation only in one direction: inhomogeneities in the cloud may lead to 
deeper penetration of the external field in a given line of sight than expected for the nominal 
value of Ay (see Andersson & Potter 2007 for detailed discussion). 

We noted in §2.3 that the p Oph cloud (Fig. 4) is subject to a much stronger ambient 
radiation field compared with Taurus, arising from the proximity of OB stars, yet polarization 
efficiencies are remarkably consistent between the two clouds at similar values of Ay. It is 
possible, of course, that differences in (e.g.) cloud structure or magnetic field geometry (§2.1) 
conspire to reduce P\/t\ in our line of sight to p Oph, and thus negate any radiation-driven 
enhancement in alignment relative to Taurus. Moreover, the radiation field near p Oph is 
enhanced primarily in the UV and may not be so different at the near infrared wavelengths 
that drive alignment of the large grains responsible for near infrared polarization in our 
model. For the largest grains, RT alignment is predicted to approach saturation in moderate 
radiation fields (Lazarian & Hoang 2007), at which point further increases in intensity will 
have no further effect on alignment. 

5. Conclusion and a suggestion for future research 

This paper yields the following primary results: 

1. Within the limits of our sample, the general trend of declining polarization efficiency 
with optical depth in quiescent regions of dense clouds is well described by a power 
law of the form p\/r x oc [Ay]~ 6 , with b « 0.52. 

2. Polarization efficiencies observed toward YSOs indicate that grain alignment tends to 
be enhanced in molecular clouds that contain regions of active massive star formation. 

3. Although our results are based on a dataset dominated by broadband polarimetry, 
observations of ice-feature polarizations in a small subset of the targets are particularly 
powerful in elucidating the efficiency of grain alignment in molecular clouds. 

4. Of extant grain alignment models, the radiative torques mechanism seems most effec- 
tive in producing degrees of alignment adequate to explain the observations. 

We note in conclusion that spectropolarimetry of ice features may provide not only 
the most stringent method of testing alignment mechanisms but also the most promising 
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technique for mapping magnetic fields internal to dense molecular clouds. It was pointed 
out some years ago that the observed reduction in polarization efficiency with extinction 
compromises the value of near infrared broadband polarimetry for mapping (e.g. Goodman 
et al. 1992, 1995). This is because the net polarization in a given line of sight tends to be 
dominated by dust in the diffuse outer layers of a cloud, where alignment is most efficient 
(the skin-depth scenario; §3.1): there is no way to separate the effect of the "skin" to reveal 
the internal field because all aligned grains in the line of sight contribute to the observed 
continuum polarization. However, the excess polarization observed in an ice absorption 
feature naturally allows this separation, because only grains internal to the cloud (where 
temperatures are low enough for ice mantles to from) will contribute to the excess. The 
amplitude and position angle of the feature component of the polarization thus traces the 
internal field. To date, the number of relevant observations is extremely small, dominated by 
bright, massive YSOs (cf. Tabled]), and includes only a single field star (Elias 16; Table 2). 
It should be possible in the future to carry out narrow-band imaging polarimetry that allows 
ice-feature polarizations and optical depths to be measured simultaneously toward large 
numbers of faint background stars over extended regions of molecular clouds. 

This research utilized data products from the Two Micron All Sky Survey, a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center, funded 
by the National Aeronautics and Space Administration (NASA) and the National Science 
Foundation (NSF). Extensive use was also made of the SIMBAD database, operated at CDS, 
Strasbourg, France. A.L. and T.H. acknowledge the support of the NSF Center for Magnetic 
Self-Organization in Laboratory and Astrophysical Plasmas and of NSF grant AST 0507164. 
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Table 1. Broadband polarization and extinction data for Taurus stars. 



Star" 


PK 

(%) 


Ap K 


Ref. 6 


A v 
(mag) 


Method 


Pk/tk 


Field stars: 
















Elias 3 


0.65 


0.02 


1. 2 


10.0 


Ej-K 


0.8 


Elias 6 


0.90 


0.08 


3 


6.8 




Ej-K 


1.6 


Elias 13 


0.85 


0.10 


3 


11.7 


Ej-K 


0.9 


Elias 15 


2.08 


0.20 


3, 4 


15.3 


Ej-K 


1.6 


Elias 16 


2.47 


0.08 


4, 5 


24.1 


Ej-K 


1.2 


Elias 19 


1.20 


0.03 


1 


2.5 




Ej-K 


5.7 


Elias 29 


0.94 


0.09 


1 


4.2 




Ej-K 


2.7 


Elias 30 


1.14 


0.04 


2 


4.9 




Ej-K 


2.8 


GJL 41919.7+263926 


0.76 


0.07 


6 


5.5: 


Ej-K 


1.6 


GJL 42003.6+264232 


0.38 


0.06 


6 


3.3: 


Ej-K 


1.4 


GJL 42014.0+263419 


0.91 


0.13 


6 


11.6: 


Ej-K 


0.9 


GJL 42041.6+263553 


0.87 


0.12 


6 


8.8 




Ej-K 


1.2 


GJL 42046.0+263313 


0.80 


0.16 


6 


6.7 




Ej-K 


1.4 


GJL 42047.6+264121 


0.51 


0.06 


6 


3.3 




Ej-K 


1.8 


GJL 42052.4+262944 


0.98 


0.12 


6 


5.5 




Ej-K 


2.1 


GJL 42055.3+264311 


0.41 


0.05 


6 


2.1 




Ej-K 


2.3 


GJL 42124.9+263323 


0.65 


0.12 


6 


4.3 




Ej-K 


1.8 


GJL 42127.8+262153 


1.29 


0.09 


6 


4.2 




Ej-K 


3.7 


GJL 42129.3+262645 


0.65 


0.08 


6 


3.0 




Ej-K 


2.6 


GJL 42134.8+262000 


1.33 


0.11 


6 


4.0 




Ej-K 


4.0 


GJL 42146.3+262205 


0.99 


0.05 


6 


3.8 




Ej-K 


3.1 


GJL 42211.6+262225 


1.01 


0.05 


6 


4.0 




Ej-K 


3.0 


GJL 42232.2+261258 


0.65 


0.12 


6 


2.7 




Ej-K 


2.9 


GJL 42234.1+262246 


1.47 


0.05 


6 


4.9 




Ej-K 


3.6 


HD 28170 


0.38 


0.07 


1, 7 


1.32 


Ev-K 


3.4 


HD 28225 


0.34 


0.02 


2 


1.22 


Ev-K 


3.3 


HD 28975 


0.82 


0.11 


1 


1.75 


Ev-K 


5.6 


HD 29333 


1.03 


0.07 


1 


1.98 


Ev-K 


6.2 


HD 29647 


0.64 


0.06 


1 


3.63 


Ev-K 


2.1 


HD 29835 


0.80 


0.06 


1 


1.14 


Ev-K 


8.4 


HD 30168 


0.59 


0.07 


1 


1.18 


Ev-K 


6.0 


HD 30675 


0.64 


0.11 


1 


1.54 


Ev-K 


5.0 


HD 279652 


0.20 


0.04 


1, 7 


0.99 


Ev-K 


2.4 


HD 279658 


0.38 


0.05 


1, 7 


1.87 


Ev-K 


2.4 


HD 283367 


0.36 


0.04 


2 


1.95 


Ev-K 


2.2 


HD 283637 


0.54 


0.05 


1 


2.28 


Ev-K 


2.8 


HD 283642 


0.42 


0.03 


2 


2.19 


Ev-K 


2.3 


HD 283643 


0.29 


0.06 


2, 7 


1.66 


Ev-K 


2.1 


HD 283701 


0.73 


0.07 


1 


2.53 


Ev-K 


3.4 


HD 283725 


1.12 


0.08 


1 


1.75 


Ev-K 


7.6 


HD 283757 


0.67 


0.09 


2 


1.70 


Ev-K 


4.7 


HD 283800 


0.69 


0.07 


1 


1.64 


Ev-K 


5.0 


HD 283809 


1.48 


0.09 


1 


5.70 


Ev-K 


3.1 
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Table 1 — Continued 



Star" 


PK 


Ap K 


Rcf. 6 


A v 
(mag) 


Method 


Pk/tr 


HD 283812 


1.11 


0.06 


1 


1.91 


Ev-K 


6.9 


HD 283815 


0.57 


0.05 


2 


1.91 


Ev-K 


3.6 


HD 283877 


0.37 


0.08 


2, 7 


0.72 


Ev-K 


6.1 


HD 283879 


0.86 


0.05 


2 


3.33 


Ev-K 


3.1 


TNS 8 


2.7 


0.3 


4 


21.5 


Ej-K 


1.5 


TNS 12 


1.5 


0.2 


4 


6.9 




Ej-K 


2.6 


TNS 13 


0.8 


0.2 


4 


2.0 




Ej-K 


4.8 


TNS 15 


1.3 


0.3 


4 


6.2 




Ej-K 


2.5 


TNS 16 


2.0 


0.2 


4 


5.3 




Ej-K 


4.5 


TNS 17 


1.2 


0.1 


4 


4.0 




Ej-K 


3.6 


Embedded YSOs: 
















Elias 1 (V892 Tau) 


1.83 


0.10 


1, 3 


13.* 


i 


Ej-K 


1.6 


Elias 5 (DG Tau) 


2.10 


0.12 


3 


5.2: 


Ej-K 


4.8 


Elias 12 (V807 Tau) 


0.90 


0.18 


3 


1.5 




Ej-K 


7.2 


Elias 18 (IC 2087 IR) 


1.48 


0.10 


3 


22.2 


Ej-K 


0.8 


GJL 41857.5+265029 


1.54 


0.11 


6 


8.4: 


Ej-K 


2.2 


GJL 41943.6+263854 c 


1.39 


0.09 


6 


7.5 




Ej-K 


2.2 


HN Tau 


2.74 


0.57 


3 


6.5: 


Ej-K 


5.0 



a Elias, GJL, and TNS identifications are from Elias 1978b, Goodman ct al. 1992, 
and Tamura et al. 1987, respectively; the Goodman et al. numbering is based on 
1950 coordinates. 

b Sourccs of polarimetry: (1) Whittct et al. 1992 with adjustment to the standard 
K passband (see Gcrakincs ct al. 1995); (2) Whittct et al. 2001; (3) Tamura & 
Sato 1989; (4) Tamura ct al. 1987; (5) Hough ct al. 1988; (6) Goodman et al. 1992; 
(7) Extrapolation based on fit to p(A). 

C IRAS 04196+2638; YSO status determined by Luhman ct al. 2006. 



Table 2. Comparison of broadband and feature data for Taurus field star Elias 16. 



Description 


A 


PA 




Px/rx 


Ref. a 




(p,m) 










K band 


2.2 


2.47 ± 0.08 


2.02 + 0.06 


1.22 + 0.08 


4, 5, 8 


H2O ice feature 


3.0 


1.26 ± 0.10 


1.43 + 0.15 


0.88 + 0.16 


5, 9 


CO ice feature 


4.67 


1.20 + 0.10 


1.06 + 0.02 


1.13 + 0.12 


10 



a References: (4)-(5) as Table 1; (8) i-T-band optical depth from the present 
work; (9) H2 0-ice optical depth averaged from Whittet ct al. 1988 and Smith 
et al. 1993; (10) Hough et al. 2008. 
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Table 3. Broadband polarization and extinction data for p Oph stars. 



Star" 


PK 


Ap K 


Ref. 6 


A v 


Method 


Pk/tr 




(%) 






(mag) 








Field stars: 
















HD 145502 Sco) 


0.26 


0.01 


13 


1.06 


Ey _ 


-K 


2.9 


HD 147084 (o Sco) 


0.91 


0.04 


12 


2.70 


Ey _ 


-K 


4.0 


HD 147283 


0.41 


0.04 


1 


2.50 


Ey. 


-K 


2.0 


HD 147888 (p Oph D) 


0.71 


0.05 


13 


1.85 


Ey _ 


-K 


4.6 


HD 147889 


1.12 


0.02 


11, 12 


4.60 


E v - 


-K 


2.9 


HD 147932 (p Oph C) 


0.87 


0.05 


13 


2.10 


Ey. 


-K 


4.9 


HD 147933/4 (p Oph A/B) 


0.58 


0.03 


13 


1.79 


Ey_ 


-K 


3.9 


HD 149757 « Oph) 


0.29 


0.02 


12 


1.06 


Ey- 


-K 


3.3 


Oph 7 (S3) 


1.30 


0.06 


11 


7.9 


Ej- 


K 


1.5 


Oph 8 (S4) 


1.88 


0.05 


11 


8.5 


Ej- 


K 


2.0 


Oph 11 (GSS15) 


1.04 


0.04 


11 


9.7 


Ej- 


K 


1.0 


Oph 16 (SR3) 


1.40 


0.08 


11, 12 


6.5 


Ej- 


K 


2.0 


Oph 25 (SI) 


2.17 


0.04 


11, 12 


13.1 


Ej- 


K 


1.5 


Oph 35 (VSSG13) 


1.65 


0.04 


11 


11.6 


E.j- 


K 


1.3 


Embedded YSOs: 
















GSS 26 


5.3 


0.4 


11 


25: 


Ej- 


K 


2.0 


HD 150193 (V2307 Oph) 


1.44 


0.03 


1 


5.9 


E.j- 


K 


2.2 


Oph 14 (V2246 Oph) 


1.85 


0.05 


14 


9.8 


Ej- 


K 


1.8 


Oph 20 (VSSG1) 


2.33 


0.12 


11 


18: 


Ej- 


K 


1.2 


Oph 22 (DoAr 24E) 


0.85 


0.04 


11 


9.2 


E.j- 


K 


0.9 


Oph 23 (S2) 


3.1 


0.2 


11 


14: 


Ej- 


K 


2.0 


Oph 24 (IRAS 16233-2409) 


2.04 


0.09 


11 


11.9 


Ej- 


K 


1.6 


Oph 27 (GSS39) 


3.6 


0.3 


11 


17: 


Ej- 


K 


2.0 


Oph 29 


6.69 


0.10 


15 


30: 


t(s 


1) 


2.7 


Oph 33 (VSSG17) 


4.2 


0.3 


11 


26: 


Ej- 


K 


1.5 


Oph 36 (VSSG14) 


2.02 


0.14 


11 


8.5: 


Ej- 


K 


2.2 



a Oph designations from Elias 1978a; GSS designations from Grasdalen, Strom & Strom 
1973; VSSG designations from Vrba et al. 1975. 

b References: (1) as Table 1; (11) Wilking ct al. 1979; (12) Wilking et al. 1980; (13) 
Wilking ct al. 1982; (14) Martin ct al. 1992; (15) Hough ct al. 1989. 
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Table 4. Summary of comparison data. 



Region / sample 


Altars 


A 


Ay 


Reference 






w 


method 




L1755 field stars 


22 


2.2 


Ej-K 


Goodman et al. 1995 


L1641 YSOs 


15 


2.2 


E>H-K 


Casali 1995 a 


H2 0-ice features in YSOs 


8 


3.0 


r(sil) 


Holloway et al. 2002 6 


CO-ice feature in YSO W33A 


1 


4.7 


T(sil) 


Chrysostomou et al. 1996 


Silicate features in YSOs 


19 


9.7 


r(sil) 


Smith et al. 2000 c 



a Only point-like sources with pn < pn were selected. 

b Cited reference includes a compilation of data from earlier literature. 

c Only sources with silicate spectra dominated by absorption were selected (see Table 0. 



Table 5. Broadband and feature polarization efficiencies compared for YSOs. a 



YSO 


A v 
(mag) 


Pk/tk 


Px/rx b 
ices 


Px/rx 
silicate 


AFGL 490 


31 


2.0 


5.1 


2.2 


AFGL 896 


68 




3.5 


1.6 


AFGL 961 


36 






1.0 


AFGL 989 


29 


1.7 


0.8 


1.6 


AFGL 2136 


77 


2.8 


1.8 


1.7 


AFGL 2591 


50 


2.6 




2.0 


AFGL 2789 


13 






2.3 


AFGL 2884 


70 


2.6 




0.8 


AFGL 4176 


72 






0.8 


M8E IR 


22 






0.8 


Mon R2 IRS2 


70 


2.2 


2.0 


1.9 


NGC 7538 IRS1 


117 


1.1 




1.0 


OMC-1 BN 


58 


3.8 


3.2 


3.9 


OMC-1 IRc4 


65 






1.2 


Oph 29 


30 


2.9 


0.8 


2.1 


SVS 13 


23 




3.2 


3.8 


W3 IRS5 


117 


1.4 




1.0 


W33A 


120 


1.5 


1.0 


1.6 


W51 IRS2 


22 






5.0 



a Values significantly enhanced relative to the trend for Tau- 
rus field stars are in bold. 

b Values of px/^x for ices refer to data for the A = 4.67 /im 
CO feature in W33A, and to the A = 3.0 fim H2O feature in 
all other objects. 
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Table 6. Key to plotting symbols used in the figures. 



Symbol 


Meaning 


Circle 


Photometric colors (Fig. 1 only) 


Square 


Broadband K (2.2 ftm) a 


Upright triangle 


H 2 0-icc (3.0 fim) 


Inverted triangle 


CO-icc (4.67 fan) 


Diamond 


Silicate (9.7 /im) 


Filled symbol (any) 


Background field star 


Open symbol (any) 


Embedded YSO 



a AlsO Amax, Pmax in Figs. 9-10. 
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Fig. 1. — JHK color-color diagram for Taurus stars (Table 1) with extinction estimates 
based on infrared data (Ej_ K method). Filled circles denote field stars, open circles denote 
YSOs. Photometry is from the Two-Micron All-Sky Survey (2MASS). Solid curves near the 
origin are intrinsic colors for normal unreddened stars of luminosity classes III and V (upper 
and lower branches). The arrow is a sample reddening vector with slope Ej_ H /E H _ K = 1.75 
(Whittet et al. 2007) and amplitude equivalent to A v ~ 5 mag. The dotted diagonal lines 
delimit the expected loci of stars with normal intrinsic colors, subject to reddening. The 
YSOs are subject to additional displacement resulting from circumstellar infrared emission 
(e.g. Itoh et al. 1996), but nevertheless show a general trend roughly parallel to the reddening 
vector. 
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Fig. 2. — Plots of polarization efficiency (p\/r\) vs. visual extinction (Ay) for Taurus stars. 
Frames (a) through (d) show 2.2/im broadband data for field stars (Table f): (a) the entire 
sample; (b) subset with Ay from the Ey-K method; (c) subset with Ay from the Ej^k 
method (no spectral types); (d) subset with Ay from the Ej_k method (with spectral 
types). Frame (e) plots 2.2/xm broadband data for Taurus YSOs (Table 1). Frame (f) plots 
values based on 3.0 /xm (H 2 0) and 4.67 /jm (CO) ice-feature spectropolarimetry (Table 2) 
for the Taurus field star Elias 16. The thick diagonal line in each frame is the least-squares 
power-law fit to the 2.2 /im data for all field stars (as in frame a), and the thinner diagonal 
lines represent the approximate upper and lower bounds of the distribution (±0.34 dex 
relative to the power-law fit). The vertical dotted line in each frame denotes the "threshold 
extinction" for H 2 0-ice in the cloud (A v = 3.3 mag; Whittet et al. 2001), above which the 
grains typically acquire mantles. 
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Fig. 3. — Taurus data compared with a simple "skin-depth" model for the variation of po- 
larization efficiency with extinction. All data plotted in Fig. 2 are included here in a single 
frame; the plotting symbols are explained in Table El The thin diagonal line is the fit to 
2.2 (Urn field-star data, as in Fig. 2, and the vertical dotted line denotes the ice threshold 
extinction. The thick diagonal line assumes that p(A)/r(A) oc [Ay] -1 , as expected if po- 
larization is effectively constant, variations in p(A)/r(A) with Ay resulting merely from the 
increase in optical depth. 
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Fig. 4. — Polarization efficiency vs. visual extinction for the p Oph dark cloud. The plotting 
symbols are explained in Table EJ i^-band data are from Table EJ H^O-ice and silicate data 
for YSO Oph 29 (Ay ~ 30 mag) are from sources listed in Table HI The diagonal lines are 
identical to those in Fig. 2, i.e., the power-law fit and upper/lower bounds to the distribution 
of field-star data for the Taurus cloud. 
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Fig. 5. — Polarization efficiency vs. visual extinction for field stars observed through the 
filamentary dark cloud L1755 (see Table H]). The solid diagonal lines are identical to those in 
Fig. 2, i.e., the power-law fit and upper/lower bounds to the distribution of field-star data 
for the Taurus cloud. The dotted line is a p(A)/r(A) oc [Ay] -1 model, as in Fig. 3. 
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Fig. 6. — Polarization efficiency vs. visual extinction for YSOs in the L1641 dark cloud 
(see Table HJ). Diagonal lines represent the power-law fit and upper /lower bounds to the 
distribution of data for Taurus field stars, as in Fig. 2. 
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Fig. 7. — Polarization efficiency vs. visual extinction for solid-state features observed in 
YSOs (see Tables H] and [5]) : H 2 and CO ice data are plotted in frame (a), silicate data in 
frame (b); i^-band continuum data as available for the same stars are plotted in frame (c). 
Diagonal lines represent the power-law fit and upper /lower bounds to the distribution of 
data for Taurus field stars, as in Fig. 2. 
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Fig. 8. — Summary plots of all data in Figs. 2-7 (field stars: frame a; YSOs: frame b). The 
plotting symbols are explained in Table [61 and the diagonal lines are again identical to those 
in Fig. 2. 
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Fig. 9. — Plot of A max vs. Ay, comparing observational data for Taurus field stars (Whittet 
et al. 2001) with a prediction for the radiative torques model described in §4 (curve). The 
Amax value is effectively a measure of the mean size of the aligned grains. 
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Fig. 10. — Polarization efficiency vs. visual extinction, comparing data for Taurus field stars 
with predictions of the radiative torques model (dotted curves, §4). Frame (a): Pk/ t k vs. 
Ay (data from Table 1); the diagonal line is the power- law fit to the data, as in Fig. 2. 
Frame (b): p ma x/Ay vs. Ay (data from Whittet et al. 2001). The curves labeled [1], [2] and 
[3] in each frame are for R* values (§4) of 0.7, 0.4 and 0.2, respectively. 



